Correlation between stress and band shift on late transition metal (111) surfaces 
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We use ab initio local stress calculations to investigate layer-by-layer ab initio stress inside late 
transition metal (111) surfaces, focusing on the origin of stress on the surface top layer. It is found 
that the band shift on each surface layer is strongly correlated with the in-plane stress. For the top 
layer, this correlation can be explained by the Friedel model. The reduction of the local d-band 
width due to the coordination reduction is the main origin of both the d-band center shift and in- 
plane tensile stress. The changes in the directional d-d bonding character analyzed by the in-plane 
and out-of-plane projected densities of states should be an additional origin of the excess tensile 
stress, except for Ag explained mainly by the Friedel model. 
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PACS numbers: 68.35.Gy, 71.15.-m, 68.47.De 

The electronic structure in a free surface is different 
from that in the bulk crystal. This is mainly due to the 
reduced coordination number (CN) of the surface atoms. 
This difference gives rise to the intrinsic stress on the 
surface because the bonding character between the sur- 
face atoms is different from that between the bulk atoms. 
Considering that stress can be defined as an energetic re- 
sponse against strain [![, we can say that the surface 
under stress is strained even if the surface bond has the 
same bond length as that in the bulk. 

Stress on metal surfaces is an important research area 
in recent surface physics pj-ticj. In particular, in-plane 
stress and strain on 4d and 5d late transition metal sur- 
faces have aroused sizable interest because of their ef- 
fects on chemisorption properties in catalysis Hi, Il2 |. 
Mavrikakis et al. [l3( showed that increased tensile stress 
leads to stronger adsorption of CO or O on Pt surfaces 
due to a shift in the center of d bands with strain. A 
bimetallic surface, i.e., a metal surface with an adlayer of 
a different metal species, is regarded as a practical instru- 
ment for producing such a strained surface fill . 

EH. The 

difference in lattice constants between the adlayer and 
substrate metals results in stress on the adlayer and its 
interfaces. Because stress in atomistic scales is not easily 
attainable in experimental methods [J [l^ , ab initio sim- 
ulation based on density-functional theory is expected to 
play an irreplaceable role for investigating the stress on 
the surface. 

In spite of the importance of stress on late transition 
metal surfaces, the origin of tensile stress on these sur- 
faces is still controversial even in ab initio theoretical 
arguments: in Ibach's model the charge redistribu- 
tion from a dangling bond to a surface bond strengthens 
the latter. Kollar et al. 0] attributed the origin to the 
depletion of sp electrons at the surface layer while d elec- 



trons are unchanged there. Fiorentini et al. 0| proposed 
a model where the depletion of antibonding d electrons 
strengthens the surface bond. To our knowledge, these 
models still coexist when explaining the origin of surface 
stress on late transition metal surfaces 0, E3] ■ 

In this paper, with a view to resolving this contradic- 
tion, we investigate the local stress distribution inside 
Ad and 5d late transition metal surfaces via the ab ini- 
tio local stress analysis developed in our preceding paper 
[TBI ] , which can reveal more details about the stress state 
on each surface layer. We then find the correlation be- 
tween the band shift and the in-plane tensile stress on 
each surface layer. We find that this correlation can be 
explained by a simple Friedel model based on the tight- 
binding second moment theory [ljj ]: for the surface top 
layers of various metals, the in-plane tensile stress can be 
explained systematically by the reduction in the d-band 
width, inducing the d-band shift. Furthermore, we show 
that the changes in the directional d-d bonding character 
observed in the in-plane and out-of-plane projected den- 
sities of states (PDOSs) can explain an additional origin 
of the excess tensile stress as discussed by Fiorentini et 
al 0. 

The calculation was performed using the projector aug- 
mented wave (PAW) method 0, El with our in-house 
ab initio code QMAS (Quantum MAterials Simulator) 
[l9l | . We use the generalized gradient approximation [2(| 
for the exchange-correlation functional. The Gaussian 
smearing method [2lT | with a smearing width of 0.1 eV 
is employed to introduce partial occupancy. This study 
deals with an fee (111) surface which has the least surface 
energy among the fee surfaces because of its high CN. 
The (lxl) surface is simulated by a slab model where 
each 14-layer slab is separated by a vacuum of six-layer 
thickness. The k-point mesh in the full Brillouin zone 
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TABLE I: Theoretical surface properties of 4d and 5d late 
transition metals: an denotes the lattice constant; and 7 and 



r are 


the surface energy 


and stress, respectively. 






an [Al 
U L J 


•7 fJ/m 2 l 


t fJ/m 2 l 

' l' J / L±1 i 


Rh 


3.848, 3.85 a 


2.02, 2.01 a 


2.61, 2.73 a 


Pd 


3.951, 3.96 a 


1.35, 1.33 a 


2.71, 2.57 a 


Ag 


4.145, 4.16° 


0.73, 0.76 a 


0.82, 0.79 a 


Ir 


3.883, 3.877 b 


2.31, 2.06 6 


4.72, 4.37 6 


Ft 


3.983, 3.978 b 


1.49, 1.49 6 


4.47, 4.25 6 


An 


4.164, 4.174 6 


0.70, 0.70 b 


2.50, 1.76 b 



a PAW, Ref. 
b PAW, Ref. 
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FIG. 1: In-plane layer-by-layer stress and energy distribution 
for a 14-layer slab of the a transition metal (111) surface. The 
data from the surface layer to the 8 th layer are plotted. Lines 
are a guide for the eyes. The in-plane stress state on an fee 
(111) surface can be expressed by a single isotropic stress due 
to the hexagonal symmetry of the surface. 



was 32x32x2 for the supcrcclls. The cutoff energy for 
the valence wave function was set to 544 eV in all the 
calculations. 

The local stress is obtained by integrating the stress 
density [22[ expressed in the PAW framework 15 1 over 
partial volumes. We can define each local region con- 
taining a single atom by setting the boundaries between 
atoms so that the gauge-dependent term originating from 
the kinetic energy density is integrated to zero in each lo- 
cal region. Then we obtain the local stress as the layer- 
by-layer stress in an fee (111) surface slab as performed in 
an Al surface fl5| . In the following, we express the layer- 
by-layer stress in unit of eV/atom. The sign of stress is 
positive for tensile stress. 

First, we discuss the layer-by-layer stress distribution 
shown in FigJT] The surface lattice constants were deter- 
mined so that the layer-by-layer stress is nearly zero at 
the center of the slab. The conventional surface energy 
or stress / is expressed as follows: 
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where g is the laycr-by-laycr energy or stress in FigJTJ S 
is the unit area of the surface, and k is the layer index 



FIG. 2: Correlation between the layer-by-layer in-plane stress 
and the shift of average band energies on relaxed and unre- 
laxed surfaces. The data from the surface layer to the 7 th 
layer are plotted. Least squares lines fitted to the data are 
also drawn. The labeled data are those on the surface top 
layers. 



as shown in Fig|TJ As shown in Table HI the obtained 
surface properties are in good agreement with previous 
PAW calculations. 

The surface stresses on late transition metal (111) sur- 
faces are tensile. However, the stress distribution is not 
limited to the surface top layer but spread to the layers 
inside the surface as shown in FigJT] For example, in the 
Pt case, the third and fourth layers from the surface top 
layer are subjected to compressive stress. This indicates 
that the conventional surface stress does not directly re- 
flect the local state of each atomic layer since it can be 
given as the sum of the layer- by- layer stress. Here we 
consider the layer- by- layer stress as a quantity express- 
ing the local state on each atomic layer at the surface. 

We next compare the layer-by-laycr stress with 
band shifts which are considered to directly affect the 
chemisorption properties of late transition metal surfaces. 
The band shift of each atomic layer is defined as the dif- 
ference in the average band energy between the surface 
layer and the bulk-like layer at the center of the slab. 
We obtain the average band energy of each atomic layer 
by integrating its local density of states (LDOS). Figure 
[H shows the results. There is a correlation between the 
layer-by-layer stress and the band shift on each atomic 
layer. For the top layers of all the late transition metal 
surfaces, we observe a clear correlation between the ten- 
sile stress and the upward band shift. This suggests that 
the origin of the stress on the surface layer can be at- 
tributed to the physics inducing the band shift. This re- 
sult coincides with the fact that the tensile strain pushes 
up the average energy of d-bands on the late transition 
metal surface Furthermore, this figure shows that 

the local stress is perturbed by the surface relaxation. 
The effect of the surface relaxation on the surface stress 
is pointed out elsewhere 0, Q . 

The band shift on each late transition metal surface is 
mainly attributed to the reduction in d-band width [23| . 
because of a constant filling of the LDOS of the surface 
top layer to keep the charge neutrality. Thus we examine 
the applicability of the Friedcl model, where the surface 
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FIG. 3: Comparison of the ab initio stress on the surface top 
layer and the stress obtained with the Friedel model. The ab 
initio stress and SX are obtained with ab initio calculation for 
the unrelaxed surface. A function a<5A 2 + fi is fitted to the 
Friedel stress data, where a and /3 are fitting parameters. 



energy of a transition metal is expressed mainly by the 
energy increase in the surface atom caused by the reduc- 
tion in the d-band width due to the CN reduction within 
the rectangular d-band model based on the tight-binding 
second moment approximation (lo| . The in-plane sur- 
face stress can be derived from the surface energy in the 
Friedel model without considering any charge redistri- 
bution. By using the energy derivative with respect to 
strain [lj], the in- plane stress a on the surface top layer 
is given as follows: 




(2) 



where E c is the bulk cohesive energy and Rj is the equi- 
librium distance to a nearest neighboring atom j. The 
parameters of the tight-binding model, p and q, arc given 
in Rcf. 16, and Z\> and Z s are the CN of the bulk and sur- 
face atoms, respectively. Experimental values are used 
for E c and Rj fHf- The sum in Eq. (J2j) runs over six 
neighboring surface atoms on an fee (111) surface. For 
convenience's sake, here we call this stress "the Friedel 
stress" . 

Figure [3] compares the Friedel stress and the ab ini- 
tio local stress. The magnitude of the Friedel stress is 
comparable to that of the ab initio stress. This indicates 
that the Friedel model can explain the main origin of the 
in-plane stress on the late transition metal surface. The 
agreement between the Friedel stress and the ab initio 
one is impressive for Ag, implying that Ag is a typical 
metal whose surface properties are well described by the 
Friedel model. 

We can assume that within the Friedel model E c in- 



creases quadratically with the bulk d-band width [25 1. 
This model also shows that the band shift, the square 
root of the second moment of d-band PDOS A, and its 
difference between the bulk and the surface SX = A s — Ab 
are all proportional to the bulk d-band width |l6} . Be- 



cause the Friedel stress is governed by E c in Eq.([2|), it 
is expected that the stress increases quadratically with 
the bulk d-band width as E c does. In FigEl a quadratic 
function of 5 A is well fitted to the Friedel stress for each 
4d or 5d surfaces because SX is proportional to the bulk 
d-band width in each metallic species. This explains why 
the ab initio stress increases with the band shift in Fig [2] 
i.e., because 5X can be regarded as the band shift. Such 
a correlation is observed only on the surface top layer 
in Fig[3] because the CN reduction occurs only on that 
layer. 

In Fig[31 however, there exist substantial differences 
between the Friedel stress and the ab initio one for the 
metals other than Ag, especially for 5d metals. This 
discrepancy is apparently caused by another mechanism 
that cannot be included in the Friedel model. One of 
the most significant effects that cannot be included in 
the Friedel model is the detailed changes in the PDOS in 
the surface atoms, other than the d-band width change. 
In particular we have to examine the changes in the di- 
rectional d-d bonding characters. Thus we have exam- 
ined the in-plane and out-of-plane PDOSs of d bands 
(d|| and d±) for the surface top layers of all the met- 
als, which were obtained as the projection to d xy and 
d x 2_ y 2, and to d zx , d yz , and d 3z 2_ r 2, respectively (The 
x, y and z axes are parallel to the crystallographic direc- 
tions (110), (112) and (111), respectively). Results are 
shown in Figd^i. We found that near the Fermi level the 
surface density of states is reduced in dy , while it is in- 
creased in dj_ compared to the bulk: the charge in the 
antibonding states is depleted from the in-plane surface 
bonds while it is not for the dangling bonds. In observing 
the shape of the entire PDOS in FigdjD [26[, we did not 
find any such depletion in either the bonding or nonbond- 
ing states. This directional d electron depletion can cause 
the surface bonds to contract and result in in-plane ten- 
sile stress on the surface atoms. Furthermore, the larger 
depletion in Ir or Pt than in Rh or Pd (FigHJi) can more 
significantly affect the tensile stress in Ir or Pt as shown 
in Figj3] Such d-electron depletion is never observed for 
Ag because its d band lies deeply below the Fermi energy. 
This is considered to be the reason the ab initio stress 
and the Friedel one show good agreement in Ag as shown 
in Fig[3] In light of the above discussion, we believe that 
the change in the directional d-d bonding characters fur- 
ther contributes to the stress on the surface top layer on 
top of the contribution described by the Friedel model. 

The d-electron depletion as the origin of surface stress 
has been discussed in several studies. In late transition 
metal surfaces, the d band shifts upward so as to keep 
the charge neutrality of the d band as mentioned above. 
Fiorentini et al. postulated that this upward shift pro- 
duces depletion of d electrons from the antibonding states 
[2| . Blanco-rey et al. argued that the excess tensile stress 
comes from a decreased compressive kinetic stress com- 
ponent induced by directional charge redistribution from 
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FIG. 4: Surface PDOSs of late transition metals: (a) dif- 
ference in the surface dii- and dx-band PDOS between the 
relaxed surface and the bulk in the vicinity of the Fermi en- 
ergy, (b) whole d|| PDOS of the surface top layer and the bulk. 
The shaded area indicates the depleted d electrons near the 
Fermi energy. 



surface bonds to dangling bonds. The PDOS results we 
obtained regarding the antibonding d-electron depletion 
in the in-plane bonding on the surface top layer are con- 
sistent with these arguments. 

The results we obtained for Au, however, do not seem 
to be simply explained by the above discussion: the ab 
initio stress is substantially larger than the Friedel one, 
while no clear antibonding d-electron depiction is ob- 
served in du . This discrepancy could come from the 
strong s-d hybridization, which governs the structure ten- 
dency of Au clusters with low CN [27, 28 1. Furthermore, 
complicated charge redistribution in Au surfaces between 
s and d surface states has been discussed [2!| . The sod 
redistribution can contribute to the stress on the surface 
top layer, especially in Au. Such rclativistic effects in 5d 
metals are not directly included in the above discussion. 



In conclusion, we revealed that the in-plane tensile 
stress and the upward d-band shift at the top layer of 
late transition metal surfaces are inseparable phenom- 
ena, because both are caused by the d-band width re- 
duction due to the CN reduction, as clearly described 
by the Friedel model. This provides us a new viewpoint 
on the chemical properties of late transition metal sur- 
faces, where the aspects of both orbital interactions and 
structural or stress relaxation should be involved in the 
chemical properties. Of course, there is an additional 
factor that contributes to the in-plane tensile stress that 



is not included in the Friedel model, i.e., the changes in 
the directional d-d bonding characters observed in the 
PDOSs of the surface top layer. The present findings 
should provide a new perspective for the properties of 
late transition metal surfaces. 
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